Dynamical arrest of ultracold lattice fermions 
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j We theoretically investigate the thermodynamics of an interacting inhomogeneous two-component 

Fermi gas in an optical lattice. Motivated by a recent experiment by L. Hackerrmiller et al, Sci- 
ence, 327, 1621 (2010), we study the effect of the interplay between thermodynamics and strong 
£Nj ' correlations on the size of the fermionic cloud. We use dynamical mean-field theory to compute the 

» ! , cloud size, which in the experiment shows an anomalous expansion behavior upon increasing attrac- 

tive interaction. We confirm this qualitative effect but, assuming adiabaticity, we find quantitative 
agreement only for weak interactions. For strong interactions we observe significant non-equilibrium 
effects which we attribute to a dynamical arrest of the particles due to increasing correlations. 

PACS numbers: 05.30.Fk, 03.75.Ss, 67.85.-d, 71.10.Fd 
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Introduction — At the heart of many condensed-matter 
phenomena lies the interplay between strong correlations 
and temperature. However, even the minimal model in- 
corporating these effects, the Hubbard model, withstands 
an exact solution. Ultracold quantum gases in optical lat- 
tices provide a new way to emulate the physics of these 
model systems in a highly controlled way not only 
can almost all system parameters be tuned with ultimate 
precision, but also the microscopic details underlying the 
model are fully known. In this way the phase diagram 
of the Bose-Hubbard model, which includes the bosonic 
Mott insulator, has been mapped out 2]. For fermionic 
lattice systems, which directly correspond to electrons in 
crystalline solid-state lattices, the formation of an incom- 
pressible Mott insulating state has been observed as well 
[H, HJ. Using Feshbach resonances, it has also been_pos- 



sible to implement the attractive Hubbard model |5j, Q , 
which sustains an s-wave superfluid state at low temper- 
ature and entropy. 

One ultimate goal of this research direction is to estab- 
lish a detailed description of the Fermi-Hubbard model 
applicable to the strongly-correlated regime of high-T c 
superconductors 0, Sj- Unfortunately, current experi- 
mental entropies are too high to observe phenomena such 
as magnetic ordering While average entropies per 
particle down to S/ksN ~ 0.5 have been demonstrated 
for fermions in pure dipole traps (lo| . in optical lattices 
only values down to S/ksN ^1 — 2 could be realized 
so far 0, [ll|, which is well above the maximum en- 
tropy w here antiferromagnetic ordering can be observed 
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However, even at the current experimental entropies, 
interesting physics emerges from the interplay between 
strong correlations and thermodynamics, as has been 
studied by several authors for repulsive interactions fl3l — 
[lij in the context of Mott insulators. Recently, Hack- 
errmiller et al. experimentally @ investigated this inter- 
play between strong correlations and entropy, focusing 



on the attractive regime. By loading a two-component 
Fermi gas into a three-dimensional (3D) optical lattice in 
the presence of a harmonic trapping potential, the size 
of the fermionic cloud was measured for different inter- 
action strengths. One would naturally expect that this 
leads to an increasing cloud size for repulsive interactions 
while an increasing attractive interaction should lead to 
a decreasing cloud size. However, a counterintuitive be- 
havior of the cloud size was observed: the cloud shrinks 
upon entering the attractive regime, but reaches a mini- 
mum at relatively small attractive interaction. For larger 
attraction the cloud expands again. 

This anomalous trend was attributed to the adiabatic 
heating effect: for sufficiently strong attractive interac- 
tion, singly occupied sites have a much higher energy 
than doubly occupied ones and become energetically ir- 
relevant. This strongly reduces the available Hilbert 
space and, at constant temperature, leads to a lower en- 
tropy. In order to keep the entropy constant, the temper- 
ature increases, leading to a higher entropy in the mo- 
tional degree of freedom and therefore to the anomalous 
expansion. 

In this Letter we theoretically investigate this ex- 
periment by means of Dynamical Mean-Field Theory 
(DMFT) simulations of the Hubbard model in the pres- 
ence of a harmonic trap. Even though our calculations 
fully include the adiabatic heating effect and also find an 
anomalous expansion of the cloud size for large attractive 
interactions, we only find an agreement between our equi- 
librium calculation and the experimental data of [6] for 
weak interactions, whereas we find a significant discrep- 
ancy for strong interactions. We attribute this to a strong 
slowing down of particle transport already for moderately 
large interactions, leading to a dynamical arrest in a state 
whose size is significantly larger than the equilibrium ra- 
dius of the particle cloud. The long-lived arrested states, 
which we find for lattice fermions, bear a str ong resem- 
blance to those of attractive colloidal glasses [17]. 
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Model and Method — We consider a two-component 
mixture of fermionic atoms, loaded into the lowest band 
of a 3D cubic optical lattice in the presence of an ex- 
ternal harmonic potential. For sufficiently deep optical 
lattices, this system is well described by the inhomoge- 
neous Fermi-Hubbard Hamiltonian 
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Here c| CT and c itT are the fermionic creation and annihi- 
lation operators and hi a — $ ia c ia is the number operator, 
where a £ {|, 1} labels the two hyperfine states. The on- 
site interaction is denoted by U, the single-atom hopping 
amplitude between nearest neighbors (i,j) is J and fi 
is the chemical potential controlling the particle num- 
ber. We consider a pancake shaped harmonic trapping 
potential of strength Vq with aspect ratio 7 and define 
rf = x\ + y\ + 7 2 zj as the squared distance of site i from 
the trap center. 

We apply DMFT to obtain the equilibrium properties 
of this Hamiltonian. Employing the assumption of a local 
self-energy, DMFT treats local quantum correlations in a 
fully non-perturbative manner [18j . The Fermi-Hubbard 
model is mapped to an effective single-site impurity prob- 
lem coupled to a non-interacting bath. Here we use ex- 
act diagonalization to obtain the self-consistent solution 
of the impurity model. To incorporate the harmonic po- 
tential, we use the local density approximation (LDA), 
which is known to reliably predict the density distribu- 
tion for fermionic systems when the trap is sufficiently 
shallow 
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2fJ . In this approximation, every site is mod- 
eled as part of a homogeneous system with local chemical 
potential ^(ri) = /i — Vb^f- We obtain the density profile 
of the fermions, from which the cloud radius is extracted. 
To make a direct comparison with the experiment, we 
employ the following definition for the cloud radius Q 



R 2 = 



1 



(i) 



where N a is the total particle number per spin state. We 
consider a balanced mixture with total particle number 
N = 2N a , such that R in Eq. (JTJ is independent of 
a. To compute the entropy per lattice site, we use the 
Maxwell relation s(/x) = J 1 ^ dfj,' d n(fx') / dT and then 
obtain the total entropy by S = J^i s (/ i ( r «))- The har- 
monic potential is characterized by the characteristic en- 
ergy E c = V rl = |Vo(37A r /87r) 2 / 3 , which is the mean 
potential energy per particle of a maximally packed state 
at the bottom of the trap with total particle number N; 
r c is the corresponding radius, which is used as the length 
scale to express the cloud size R/r c . We note that in LDA 
the rescaled radius R/r c is fully determined by the aver- 
age entropy per particle S/ksN and the ratios of tunnel- 
ing to trap strength J / E c and interaction to tunneling 
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FIG. 1. Comparison of experimental cloud radius with DMFT 
results [2J]. While we find good agreement for small inter- 
actions, the experimentally observed radii are significantly 
higher than expected from theory for stronger interactions. 
We attribute this to the dynamical arrest discussed in the 
main text. Also shown are the results of a 6th order high 
temperature series exp ansion (HTSE) of the Hubbard model 
(dashed lines) [U [2a |. Since temperature decreases with in- 
creasing J/E c and U/J (cf. Fig.[2](left)) the HTSE is only ap- 
plicable to sufficiently small or negative U/J and small J/E c . 



U/J, implying that all results shown in Fig. [T][5] are in- 
dependent of the details of the trap. The experiment has 
been performed for N a rj 1.5 x 10 5 particles at a tem- 
perature T/Tp = 0.12 ±0.03 before loading of the lattice 
(where Tp is the Fermi temperature) and an entropy per 
particle of S/ksN = 1.15±0.25. In the experiment, both 
S/kjsN and J/E c were known only with a considerable 
uncertainty [6[. Since the effect of the experimental un- 
certainty in J I E c on the cloud size is larger, we chose to 
perform all calculations for a fixed entropy per particle 
of S/ksN = 1.15 and determined the value of J/E c by 
comparing the experimental radius for U/J = with ex- 
act diagonalization (ED) results. In the non-interacting 
case non-equilibrium effects are minimal, so that the ex- 
perimental data at U/J = can be reliably fitted to 
theoretical calculations [21j . The resulting J/E c agree 
with the experimental data within errorbars for low con- 
finements, and only deviate at high confinements, where 
also some heating was observed in the experiment (cf. 
supporting online material to [f|). In addition we find 
satisfying quantitative agreement with our DMFT cal- 
culations for the fixed J/E c (without additional fitting 
parameters) in the regime of weak interactions. 

Results — In Fig. [T] we show our DMFT results for the 
rescaled cloud radius R/r c versus U/J for various ra- 
tios of tunneling to trap strength J/E c , which are com- 
pared directly to the experimental data points Q. Also 
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FIG. 2. Left: Temperature vs interaction for an entropy per 
particle of S/ksN — 1.15 and different external harmonic 
confinements. Right: Cloud radius for different entropies 
S/ksN — 0.4— 1.6. The anomalous expansion becomes more 
pronounced for higher entropies. 



a comparison with a 6th order high temperature series 
expansion is included in this figure, which agrees very 
well with the DMFT calculations. While we included 
the possibility of s-wave superfluid order in our calcula- 
tions, we found that the experimental entropy is too high 
for superfluidity to be present. 

For weak interactions \U/J\ < 2, we observe very 
good quantitative agreement with the experimental data 
points, evidencing that the DMFT solution of the single- 
band Hubbard model incorporates the relevant physics 
in this parameter regime. For stronger attractive in- 
teractions we also find the anomalous expansion of the 
cloud size in the DMFT calculations. By plotting the 
temperature at constant entropy, we can directly verify 
that this is indeed due to adiabatic heating originating 
from a reduced available phase space. This is displayed 
in Fig. [2] (left) and indeed shows a strong rise for large 
attraction. When the entropy is decreased, the anoma- 
lous expansion is expected to disappear, since essentially 
ground state properties are probed in this case. This is 
clearly visible in Fig. [5] (right), where the radius for dif- 
ferent entropies is shown. For repulsive interactions an 
adiabatic cooling effect is observed, similar to what was 
found in [l2| and Q (supporting online material). Al- 
though for strong repulsion a similar reduction in phase 
space occurs, because doubly occupied sites are energeti- 
cally unfavored, this is more than compensated by addi- 
tional (spin) configurational entropy for the fermions at 
the wings of the clouds. 

Although our simulations thus fully include the adia- 
batic heating effect, our results show that the anomalous 
expansion seen in the experiment cannot be explained by 
adiabatic heating alone: At stronger interactions we find 
significant deviations from the experimentally observed 
cloud sizes, which are much larger than theoretically ex- 
pected. By measuring the temperature after unloading 
the atoms again from the lattice it was excluded that 



the discrepancy is due to heating [6j. We attribute this 
difference to strong non-equilibrium effects in the load- 
ing: after evaporative cooling, the cloud size in the pure 
harmonic trap is significantly larger than the final equi- 
librium size in the lattice. During the ramp-up of the lat- 
tice, the effective mass of the atoms increases due to a re- 
duction of kinetic energy in the lowest band of the lattice. 
This leads to a shrinking of the cloud (see Fig. [3]) under 
adiabatic conditions. The scattering length, which is ad- 
justed by using a Feshbach resonance before the loading 
of the lattice, gives rise to an interaction strength UJ J 
which rapidly grows with the lattice depth. As recently 
demonstrated [2l| , even moderate interactions slow down 
the atoms severely, thereby prohibiting the large-scale 
particle transport necessary for following the adiabatic 
path. This can be seen e.g. by comparing the initial 
radius before the ramp-up with the radius in the lattice 
(Fig. 4 (left)): In the experiment, the radius shrinks al- 
ways by the same percentage independently of the trap 
strength, which would not be the case if the cloud stays 
in equilibrium. In contrast to the situation in 21], where 
a symmetric reduction of particle transport for attractive 
and repulsive interactions was observed, we observe much 
weaker non-adiabatic effects for the repulsive regime in 
this case. This is probably due to the fact that the differ- 
ence between the radius before the ramp-up and the equi- 
librium radius at the final lattice depth is much smaller 
in the repulsive than in the attractive regime. Once this 
effect is scaled out (see Fig. [4] (right)) we see a significant 
deviation on both sides. Moreover, repulsive interactions 
decrease the number of double occupancies while attrac- 
tive interactions increase it, thereby giving rise to differ- 
ent dynamics. Interestingly, we see that the deviation in 
Fig. [4] (right) becomes independent of U / J for sufficiently 
strong attraction. This is because the experimental and 
theoretical curves in Fig. Q] - although offset - have the 
same slope. 

To quantify the proposed scenario, i.e. that the 
disagreement between the experimental and theoretical 
cloud size has to be attributed to a dynamical arrest of 
the particles due to the slowing down of particle trans- 
port, we substitute this gradual change in adiabaticity 
by a step-wise model where up to a critical lattice depth 
Vic the system is in thermal and chemical equilibrium and 
experiences fully adiabatic changes of its state. When the 
lattice depth exceeds a critical value we assume the cloud 
to be completely static, as illustrated in Fig. [3] (right). 

Within this model we only need equilibrium simula- 
tions to determine the critical lattice depth. As seen be- 
fore, within LDA the rescaled radius R = R/r c depends 
only onS = S/k B N, J = J/E c and U = U/J. This 
allows us to determine the lattice depth at which the sys- 
tem freezes by calculating J = J(V/) and U = U(Vi,a) 
as a function of the lattice depth Vi = Vi /Er using band 
structure calculations. Here a = a/ao (arj: Bohr radius) 
denotes the scattering length. We take into account that 
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FIG. 3. Interacting fermions during the lattice ramp up: In 
the experiment, the harmonic confinement E c and the scatter- 
ing length are set in the dipole trap and remain fixed during 
the lattice ramp. Left: The tunneling and interaction dur- 
ing the ramp. Right: Equilibrium cloud size for attractive 
and repulsive interactions during the ramp. The dashed line 
denotes the non-adiabatic path of the dynamical arrest. 



E c also depends on the lattice depth due to the anti- 
confining effect of the optical lattice. The critical lattice 
depth Vic is then denned as the value where the calcu- 
lated equilibrium radius coincides with the experimen- 
tally measured one, i.e.: 



R exp = R(S,J(V lc ),U(Vi c ,a)), 



(2) 



where R cxp is the experimentally measured radius and 
the function R is calculated with DMFT. 

The resulting critical lattice depths V/ c for dynamical 
arrest of the cloud are shown in Fig. [5J First, we observe 
a strong dependence on the scattering length, which is 
indeed fully consistent with the dynamical arrest hypoth- 
esis, as the scattering length linearly affects U / J, which 
determines the transport properties. We note that for 
stronger interactions the critical lattice depth is rather 
low, and almost leaves the validity regime of the single- 
band Hubbard model. Second, we observe that Vi c is 
relatively independent of the external harmonic confine- 
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FIG. 4. Left: Final radius R of the cloud relative to the 
initial radius i?i n i before ramp-up of lattice, (symbols: ex- 
periment, lines: equilibrium theory, same colors denote same 
J/E c ) Right: Relative difference between expected (AR cq = 



Rcq — Rim) and observed (AR C 
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.Rim) change in 



cloud size during the lattice loading from the expected change 
in cloud size during the ramp. 
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FIG. 5. Critical lattice depth calculated from the experi- 
mentally measured radius via Eq. as a function of the 
scattering length. The dashed line indicates the final maxi- 
mum lattice depth at Vi / Er = 7 to which the optical lattice 
was ramped up. The values for J/E c given in the legend are 
reached for lattice depth 7Er. 



ment Vq, such that the data points for different J/E c 
collapse almost to a single curve. Indeed, while a tighter 
trap increases the density and thereby slows down the 
dynamics, it at the same time also decreases the distance 
for the required particle transport. 

Conclusion — We calculated the cloud size of an in- 
teracting Fermi gas by means of DMFT and obtained 
very good quantitative agreement with experimental data 
for weak interactions. Despite qualitative agreement, 
we observe a significant discrepancy with the exper- 
iment for strong interactions, which we attribute to 
non-equilibrium effects caused by an interaction-induced 
slowing down of particle transport. The system is there- 
fore dynamically arrested at a critical lattice depth, 
which we observe to be almost independent of the har- 
monic trap. Our findings have far-reaching consequences 
for future experiments on strongly interacting Fermi 
gases in optical lattices, for instance with the goal of ob- 
serving s-wave superfluidity or antiferromagnetic order: 
According to the scenario described above, any loading 
sequence into a deep lattice should be tailored such that 
it minimizes the required density redistribution, other- 
wise the system will always freeze in a non-equilibrium 
arrested state, rendering the observation of the desired 
low-temperature properties challenging. 
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